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' Besides the opportunity for discovering new neutrino physics, solar neu- 

(N ! 

, trino measurements provide a sensitive probe of the solar interior, and thus 

■ a rigorous test of solar model predictions. We present model independent 



determinations of the neutrino spectrum by using relevant flux components 
as free parameters subject only to the luminosity constraint. (1) Without the 



4h 

Mikheyev-Smirnov-Wolfenstein (MSW) effect, the best fit for the combined 

Q_i' data is poor. Furthermore, the data indicate a severe suppression of the 7 Be 

<D : R 

• flux relative to the B, contradicting both standard and nonstandard solar 



models in general; the pp flux takes its maximum value allowed by the lu- 
minosity constraint. This pathology consistently appears even if we ignore 
any one of the three data. (2) In the presence of the two-flavor MSW effect, 
the current constraint on the initial 8 B flux is weak, but consistent with the 
SSM and sufficient to exclude nonstandard models with small 8 B fluxes. No 
meaningful constraint is obtained for the other fluxes. In the future, even 
allowing MSW, the 8 B and 7 Be fluxes can be determined at the ±(15 — 20)% 
level, making competing solar models distinguishable. We emphasize that the 
neutral current sensitivity for 7 Be neutrinos in BOREXINO, HELLAZ, and 
HERON is essential for determining the initial fluxes. The constraints on the 
MSW parameters in the model independent analysis are also discussed. 
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I. INTRODUCTION 



The solar neutrino deficit, confirmed by all existing experiments, challenges our under- 
standing of the Sun as well as of neutrinos. The purpose of this paper is to consider the 
possibility of model independent determinations of the principle neutrino flux components 
from the solar neutrino data, both at present and in the future and with and without new 
neutrino properties. These fluxes can then be compared with the prediction of any solar 
model, standard or nonstandard. In fact, if the present experiments are correct, purely 
astrophysical explanations for the flux deficit are highly unlikely: 

• The standard solar models (SSMs) |||9| are excluded by the data as summarized in 
Table |. The discrepancy cannot be reconciled by simply changing input parameters 
in the SSM calculations |L0] . 



• The lower observed rate of Homestake relative to Kamiokande is incompatible with 
astrophysical solutions [] in general. This is a much more serious difficulty than the 
simple deficit of observed neutrinos relative to the SSM expectations. Q A model 
independent analysis |TTJ suggests a complete elimination of the 7 Be flux and, in addi- 



tion, a larger depletion of the 8 B spectrum at lower energies and/or additional neutral 
current events from or v T in Kamiokande. The larger suppression of the 7 Be than 
the 8 B flux contradicts nonstandard solar models in general, including ad hoc ones. A 
distortion of the 8 B energy spectrum cannot be caused by astrophysical effects at the 



observable level [12|. Vn and v T can interact through the neutral currents in electron 



scattering in Kamiokande, and their existence in the solar flux signifies neutrino flavor 



By astrophysical solutions, we include those involving nuclear reaction cross sections in the Sun, 
but not the chlorine and gallium detector cross sections. 

2 In fact, the discrepancy in the relative rate is aggravated in models in which the 8 B flux is reduced 
(e.g., by lowering the core temperature or reducing the 7 Be(p, 7) 8 B cross section) to explain the 
Kamiokande data. 
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oscillations. 

• The problem of the larger suppression of the 7 Be flux relative to 8 B remains even if 
we ignore any one of the three data. In this sense, the data are consistent with each 
other. In particular, if we consider the Kamiokande and the gallium results only, the 
nonstandard solar models consistent with the Kamiokande result generally predict a 
gallium rate larger than 100 SNU, f\ inconsistent with the combined result of SAGE 
and GALLEX (77 ± 9 SNU). 

With standard neutrino physics, the current situation forces us to consider a serious 
problem with two or more of the experiments and a drastic revision of the SSM calcula- 
tion unless all of the experiments are wrong. The Mikheyev-Smirnov-Wolfenstein (MSW) 
mechanism [[DJ, on the other hand, provides a complete description of the data and is also 
consistent with the SSM (see Jl4[ and references therein). Because of the consistency with 



the experiments and the simplicity of the theory, we consider the two-flavor MSW solutions 
as the most attractive scenario among many proposed particle physics solutions. 

For solar astronomy, whether new neutrino physics is present or not, the central issue is 
the determination of the solar neutrino spectrum. The theory of the Sun, which is the best 
measured main sequence star, is the keystone of our understanding of stellar structure and 
evolution. Solar neutrinos are a direct, sensitive probe of the solar core, and the neutrino 
flux measurements provide an opportunity for rigorous tests of solar models, standard or 
nonstandard. 

For the SSM, the neutrino spectrum is a diagnostic of the underlying assumptions in the 
theory. The flux prediction depends on the input physics, such as the opacity calculation 
and the nuclear cross sections, whose uncertainties might be underestimated. In particular, 
the p( 7 Be, 8 B)7 cross section, which is directly proportional to the 8 B flux, was recently 
measured using the Coulomb dissociation method in the RIKEN experiment ]T3| . Although 



3 SNU (solar neutrino unit) = 1/ 10 36 atoms /sec. 
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the measurement uncertainty is still large, the preliminary result suggests the cross section 
can be 25% lower than the current standard value ||16|| . The SSM also includes simplifi- 
cations such as the omission of rotations, magnetic fields, and the gravitational settling of 
various elements. Those effects on the neutrino flux have never been quantified in the SSM 
uncertainties. 

The nonstandard solar models, most of which are constructed to explain the solar neu- 
trino deficit, assume nonstandard input parameters or nonstandard mechanisms. Examples 
are the low central temperature (Tc), low opacity, low Z, large Sn, [] large S33, small 634, 
small Sn, mixing, and weakly interacting massive particle (WIMP) models. The neutrino 
data should test the validity of such (often ad hoc) assumptions. 

To determine the solar neutrino spectrum from the experiments, one needs to extract 
from the data the magnitude of the flux, component by component. The Kamiokande 
experiment measures the 8 B flux exclusively. The radio-chemical detectors measure the flux 
components only as a weighted sum according to the energy dependence of the detector 
cross sections: the Homestake chlorine experiment is sensitive mainly to the 8 B flux, but 
also to the 7 Be, CNO, and pep fluxes; the gallium experiments measure all components, 
including the dominant pp flux. In the future, the Sudbury Neutrino Observatory (SNO) 
17| , Super-Kamiokande and ICARUS Jl9| will measure the 8 B flux with a high precision. 



BOREXINO [^] will be capable of measuring the 7 Be line spectrum. HELL AZ [ 21[ , and 
HERON |22[ will observe the 7 Be flux and the main pp flux individually. 



This solar neutrino spectroscopy can be complicated if new particle physics effects are 
present, since those effects are often energy dependent and therefore distort the energy spec- 
trum. Uncertainties in the neutrino parameters contribute to uncertainties in the neutrino 
flux, and vice versa. In the presence of the MSW effect, for example, the determination 
of the initial (undistorted) flux components requires a knowledge of the neutrino parame- 



Sn, S33, S34, and Sn are the S factors proportional to the cross sections for p+p — > 2 H+e + + z^ e , 
3 He + 3 He -> 4 He + 2p, 3 He + 4 He 7 Be + 7, and p + 7 Be -> 8 B + 7, respectively. 
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ters and, in turn, the determination of the neutrino parameters depends on the initial flux 
magnitudes. 

To extricate the neutrino flux components from the data and distinguish various compet- 
ing solar models, it is best to consider a simple and general theoretical framework including 
all standard and nonstandard solar models. Such an analysis scheme should be viable with 
and without particle physics effects. 

Variations of solar models have usually been considered in model dependent frameworks. 



Monte Carlo SSMs |[23i|24| , |10f were obtained from various input parameters normally dis- 
tributed about their most probable values. Those solar models are, however, calculated 
within the SSM and do not address the possibility of nonstandard processes omitted in the 
modeling or the possibility of input parameters grossly different from the standard values. 
The low Tc model fl25|-f27jj parameterizes the neutrino fluxes by nonstandard core temper- 



atures as power laws | 23|j24]| . The description is more general than the Monte Carlo SSMs 



since it includes a large class of nonstandard solar models. Again, however, the Tc descrip- 
tion is model dependent: there are nonstandard solar models that cannot be parametrized 
simply by nonstandard T c , such as those with nonstandard Si 7 or S34 values. 

In this paper, we consider a model independent description of solar models, characterized 
by the magnitude of each of the neutrino flux components. By setting up an analysis scheme 
as general as possible, we depart from particular theoretical constraints. We hope that the 
experiments will distinguish standard and nonstandard models and eventually identify the 
correct solar model. The purpose of this paper is to demonstrate that such a description is 
feasible and is a powerful tool in analyzing the solar neutrino data, especially once the high 
precision data from the next generation experiments are available. 

Our model independent analysis originates in Ref. ]TT[ (see also Ref. | 28|J57fl ). Here we 



elaborate the analysis and extend it to the case in which the two-flavor MSW effect is 
present. (Of course, the analysis can be generalized in the presence of any particle physics 
effect.) We consider magnitudes of the four prominent flux components, pp, 7 Be, 8 B, and 
CNO (the sum of 13 N and 15 0), as free parameters in fitting data. In doing so, we make 



minimal assumptions: 



• The Sun is in a quasi-static state, and the solar luminosity is generated by the ordinary 
nuclear reactions of the pp and CNO chains. This imposes a relation among the fluxes: 

<f>(pp) + <p(pep) + 0.958 0(Be) + 0.955 0(CNO) = 6.57 x 10 10 cm" V 1 , (1) 

where 0(CNO) is the sum of the 13 N and 15 fluxes, which are varied with the same 
scale factor. 

• Astrophysical mechanisms do not distort the shape of the energy spectrum of the indi- 
vidual flux component at the observable level. It was shown that possible distortions of 
the spectrum due to such astrophysical effects as gravitational red-shifts and thermal 



fluctuations are completely negligible ||12|| . On the other hand, particle physics effects 
such as the MSW mechanism are in general energy dependent and lead to significant 
spectral distortions. 



The detector cross section calculations [p3|,p4[p9| are correct. 



• The minor fluxes (pep (], 17 F, and hep f\ ) are set to the SSM values. 

The simultaneous (global) analysis of all data is essential in obtaining constraints on 
fluxes. No one experiment provides enough information to determine the entire neutrino 
spectrum; only by combining various experiments with different energy thresholds is it 
possible to determine each flux component and test solar model predictions. For example, 
if there are no new particle physics effects, then by combining the 8 B flux measured in 



5 The pep neutrinos are from the reaction p + e +p — > 2 H + u e . The pep flux is the largest among 
the three minor fluxes. It is strongly correlated with the pp flux in many of the nonstandard solar 



models [30 and does not vary significantly from model to model. Of course, one can also use the 



pep flux as a free parameter. 
6 The hep neutrinos are from the reaction 3 He + p — > 4 He + e + + v e . 
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Kamiokande with the Homestake and gallium experiments, one can deduce an absence of 
the 7 Be flux and a detection of the pp flux. Allowing for MSW or other particle physics 
effects, a global analysis is even more essential because one must simultaneously determine 
the initial fluxes and the MSW-induced flux reduction and spectral distortions. 

In constraining the fluxes and testing solar models, a consistent joint analysis is impor- 
tant. Taking overlaps of parameter space allowed by different experiments does not yield 
a correct estimation of uncertainties |L4]; one needs to carry out proper joint x 2 analyses, 
which are essentially identical to the maximum likelihood method in gaussian cases. Since 
we are testing theoretical models statistically, the experimental and theoretical uncertainties 
have to be incorporated; a proper treatment of the correlations among uncertainties is also 
important |fl4|j . 

When the MSW effect is present, it is best to incorporate available energy spectrum 
and day-night asymmetry data to obtain additional constraints on the neutrino parameters 
and therefore on the fluxes. Those constraints from future high-counting experiments would 
be especially useful. In this paper, however, we do not incorporate the existing energy 
spectrum and day-night data from Kamiokande ||; their uncertainties are large and do not 
significantly change the results obtained from the averaged data. [] 

The rest of the paper is organized as the following. In Section y, the constraints on 
the fluxes without introducing particle physics effects are obtained. The joint fit of the 
combined the Homestake, Kamiokande, SAGE, and GALLEX data yields a poor \ 2 value. 
The best fit is, in fact, obtained for a negative 7 Be flux, suggesting a serious problem with 
the experiments or the existence of new particle physics effects: a distortion of the 8 B 
spectrum and/or neutral current contributions from z/ M or v T in Kamiokande. Even if we 



7 Also, Kamiokande has only presented the spectrum and angular (direction of the Sun with 
respect to the nadir) distributions separately. Since they are based on the same data, the two 
distributions are correlated and cannot be used simultaneously. It is recommended that in the 
future the data be presented in bins of definite energy and angle. 
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accept this poor fit assuming standard neutrinos, the constraint on the fluxes contradicts 
nonstandard solar models in general. The constraints on the fluxes will be displayed in the 
7 Be- 8 B, pp- 7 Be, and pp- 8 B planes. (For simplicity, we refer to neutrino fluxes in units of 
the reference fluxes listed in Table [TT] unless otherwise mentioned. Those fluxes correspond 
to the Bahcall-Pinsonneault fluxes with the helium diffusion effect ||). The results when 
one of the three experiments are omitted will be also given. By ignoring one experiment, 
the uncertainties in the flux constraints become larger, but the constraints are consistent 
with those obtained from all data, and again contradict astrophysical/nuclear solutions. 

In Section |J, we also discuss possible constraints on the fluxes from the next generation 
experiments. In fact, if no new particle effects are present, SNO and Super-Kamiokande will 
determine the initial 8 B flux with a high precision, and BOREXINO, HELLAZ, and HERON 
will measure the initial 7 Be flux exclusively. The flux constraints from hypothetical results 
from these experiments with various central values and various measurement uncertainties 
will be examined. 

In Section |TT|, we consider the constraints on the fluxes when the two-flavor MSW ef- 
fect for transitions v e — > or v e — > v T is present. Our analysis scheme can in principle 
be applied to other particle physics scenarios, such as three-flavor MSW, transitions into 
sterile neutrinos, vacuum oscillations, a large neutrino magnetic moment, neutrino decays, 
flavor changing neutral currents, violation of the equivalence principle, etc. We consider the 
two-flavor MSW solution because of its simplicity and viability. It is likely that if two-flavor 
MSW is indeed occurring, there will be enough complementary information [e.g., from spec- 
tral distortions, day-night asymmetries, and SNO neutral current (NC) measurements] to 
establish it as the most likely candidate even allowing nonstandard solar models [^4[ . Of 
course one could never rigorously exclude the possibility of more complicated scenarios, such 
as the simultaneous importance of transitions into (or u T ) and sterile neutrinos, which 
would interfere with the model independent flux determinations. For this, one must invoke 
Occam's razor. 
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Once the MSW parameters are introduced as additional free parameters in the joint fit, 
constraining the fluxes from the data is not trivial. The MSW effect can distort the energy 
spectrum depending on the parameters, and can change the contribution from different flux 
components. With the existing data, we can constrain the 8 B flux only roughly. Even 
though the chlorine and gallium experiments have a sensitivity to the 7 Be flux, the survival 
probability of the flux can be zero for the MSW small-angle (nonadiabatic) solution, and 
no meaningful constraint is obtained for the Be flux. To constrain the fluxes and the 
MSW parameters simultaneously, we need results from the future experiments, especially 
the neutral current measurement in SNO and the 7 Be neutrino measurement in BOREXINO, 
HELLAZ, or HERON. The neutrino-electron scattering mode in these 7 Be measurements has 
a sensitivity to the neutral current interactions with and u T , whose cross sections are 21% 
of z/ e 's at this energy. We will present the possible constraints assuming various outcomes 
from those experiments, and show that such a model independent analysis can determine the 
solar neutrino spectrum with an accuracy sufficient to test solar model predictions. We note 
that our choice of the hypothetical results from the SNO NC and BOREXINO experiments 
are minimal; additional information from the SNO charged current (CC) rate, the Super- 
Kamiokande rate, and the spectral and day-night asymmetry measurements in SNO and 
Super-Kamiokande should make the constraints even better. 

II. FLUX CONSTRAINTS ASSUMING STANDARD NEUTRINOS 

A. Present 

We consider the constraints on the neutrino fluxes from the updated solar neutrino data 
listed in Table D. The main results are displayed in the 7 Be- 8 B plane (although some of the 
results are also shown in the pp- 7 Be and pp- 8 B planes). When the data are fit, the the pp 
and CNO fluxes are varied freely for each 0(Be) and 0(B), subject only to the luminosity 
constraint. This representation in the 7 Be- 8 B plane is effective since it can display every 
possible solar model, standard or nonstandard, that satisfies our minimal assumptions. Since 
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predictions for those fluxes vary substantially from model to model, the 7 Be- 8 B plane also 
provides a useful diagnostic for experimentally distinguishing competing solar models. 

When the Kamiokande, Homestake, and the combined gallium experiments of SAGE 
and GALLEX are fit separately, the constraints on the 7 Be and 8 B fluxes are shown in 
Fig. [I]. The fits include the uncertainties in the radio- chemical detector cross sections and 
in the minor fluxes, which are set to the SSM values. The Kamiokande result constrains the 
8 B flux; the Homestake result constrains the Be, 8 B, and CNO fluxes; the gallium results 
constrain all fluxes including the pp. 

When all data are fit simultaneously, the allowed fluxes are severely constrained, as 
shown in Fig. |2|. The best fit for physical (i.e., non- negative) fluxes are obtained for zero 
7 Be and CNO fluxes, and the 8 B flux is about 40% of the SSM prediction; the absence of 
the 7 Be and CNO fluxes forces the pp flux to be the maximum value (1.095 SSM) allowed by 
the luminosity constraint (Eq. p]). These constraints at la uncertainties are summarized in 
Table fTJ; they are also listed as absolute fluxes in Table [TV]. This model independent result 
displays serious problems for any purely astrophysical explanation for the solar neutrino 
deficit fTJ: 

• The best fit is poor; the \ 2 minimum is in fact obtained for the unphysical value 
0(Be)/0(Be)ssM = —0.5. Imposing positivity of the flux, xLin = 3.3/1 d.f., Q which is 
excluded at 93% confidence level (C.L.) That is, any possible solar model explanation 
consistent with our minimal assumptions is excluded at least at the 93% C.L. 

• Even if one accepts this poor fit, the allowed fluxes are difficult to explain. Since 8 B 
is produced through the reaction p + 7 Be — > 8 B + 7, any reduction in 7 Be causes at 
least an equal reduction in 8 B. Therefore, unless there is some independent mechanism 

8 The fit is in fact for d.f. [3 data - (4 parameters - 1 constraint)]. For x 2 values other than 
zero, there is no standard statistical interpretation exists other than to conclude that this model is 
excluded. To quantify the confidence level, we allow 1 d.f. by considering that the 7 Be flux is fixed 
to zero. 
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to suppress only the 7 Be flux [] or the uncertainty in the 7 Be electron capture rate is 
grossly underestimated, the 8 B flux is expected to be reduced more than the 7 Be flux, 
contrary to the data. 

• Finally, various standard and nonstandard models are also displayed in Fig. 0: the 
Bahcall-Pinsonneault SSM including the helium diffusion effect H, the Bahcall-Ulrich 



1000 Monte Carlo SSMs |23|, Turck-Chieze-Lopes SSM @, the low Z model |23|,|4 



the low opacity models with the opacity reduced by 10 and 20% [31], the WIMP 



model []3^| , the large Sn models f33f , the small S34 model ||, the large S33 model 



the mixing models |34]], the Dar-Shaviv SSM |35|| , and the high Y model |23| . |24 



Also shown are models parametrized by a lower Tq (which approximately incorporates 
many of the explicit models) and a lower 5*17. As seen in Fig. |2], none of those solar 
model predictions are even close to the observations. 



We also note that a lower Sn value, suggested by the RIKEN experiment [15[], aggravates 
the problem with astrophysical/nuclear solutions, contrary to the general notion. A lower Sn 
value can make the theory prediction for the 8 B flux smaller and closer to the Kamiokande 
result, which leaves little room to introduce other astrophysical/nuclear effects (e.g., a lower 
Tc) to reduce the 7 Be flux, failing to explain either the Homestake or the gallium results. 

This complete phenomenological failure of astrophysical solutions suggests nonstandard 
particle physics effect such as the MSW effect, or serious problems with the experiments 

0- 

Even if only the Kamiokande and gallium results are considered, there is still essentially 
no viable theoretical explanation. Although the best fit somewhat improves (x 2 /0 d.f. = 
1.2), the obtained fluxes displayed in Fig. are consistent with a complete depletion of the 
7 Be flux, while the 8 B flux is about half of the SSM prediction (Table [ITT] and [TV]). This is 



9 For example, both 4>(Be) and 0(B) could be suppressed by a low Tc, and (p(B) could then be 
enhanced by a larger Sn- However, for any realistic £17, this enhancement would be negligible. 
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again in severe contradiction with nonstandard solar models in general. The nonstandard 
solar models that are significantly inside the 99% CL. contour in Fig. ^(a) are the small 
5*34 model, the large S33 model, and ad hoc mixing models that involve a core with 0.4 
and 0.8 solar masses that is mixed continuously. These models also predict non-zero CNO 
fluxes, while the CL. contours in Fig. [5] corresponds to zero CNO flux. The non-zero CNO 
contribution further aggravates the disagreement. When the CNO flux is fixed to the SSM 
value, the constraint for the combined Kamiokande and gallium results is displayed in Fig. £|. 

The discrepancy between solar model solutions and the combined Kamiokande and gal- 
lium result can be described in another way. Nonstandard models yield a wide variety of 
fluxes, and therefore a large range for their gallium predictions: from the 78 SNU of the 
luminosity limit Q to 303 SNU of the maximum rate model ||. However, the 8 B flux, which 
has the largest uncertainty among the major fluxes, has been constrained by Kamiokande, 
and this in turn constrains the gallium predictions of nonstandard models. Such a constraint 



was considered in the SSM framework with the Monte Carlo method |TIJ, but here we allow 
nonstandard models as well. Displayed in Fig. ^ along with the gallium data are the gal- 
lium predictions of various SSMs and also of nonstandard solar models that are consistent 
with or close to the 8 B flux observed in Kamiokande: the model with Su normalized to the 
Kamiokande result, the low Tc model with a reduction of 4%, the model with a larger Sn 
(30%) and a lower T c (5%), the low S u (50%) model, the Dar-Shaviv SSM @, the low 
opacity model f3Tf , the large Sn models that predict 0(B)/0(B) SSM = 0.39 and 0.57 |33 



and mixing models f3~4]| . The uncertainties include the 8 B uncertainty due to the Kamio- 



kande uncertainty (14%), but the dominant contribution is from the gallium cross section 
uncertainty. From this list, we obtain 

Gallium rate consistent with Kamiokande > 100 SNU, (2) 



10 This corresponds to zero 7 Be, CNO, pep, and 8 B fluxes and the pp flux with the maximum 
value allowed by the luminosity constraint (Eq. pi). 
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while the combined gallium rate of SAGE and GALLEX is 77 ± 9 SNU. 

The lower limit of 100 SNU can be roughly understood as followings. The nonstandard 
solar models considered here all predict smaller reductions of the 7 Be than the 8 B flux, and 
this, combined with the Kamiokande result, gives the lower limit on the 7 Be flux to be about 
half of the SSM value, contributing at least 18.3 SNU to the gallium rate. The pp and pep 
fluxes do not depend significantly on solar models; the luminosity constraint and a decrease 
in the Be flux result in an increase in the pp flux by 5% (3.7 SNU). Adding these (pp, pep, 
7 Be, and 8 B) gives a total of 102.9 SNU, with uncertainties in the treatment of the CNO 
fluxes, the gallium detector cross section, and the 8 B measurement in Kamiokande. 

This discrepancy is extremely important because it is independent of the Homestake 
result, but displays exactly the same symptom as in the Kamiokande- Homestake compari- 
son: the absence of the 7 Be flux, for which astrophysics offers no explanation. Furthermore, 
experimental developments in the near future will significantly influence the situation. The 
calibration of the gallium detectors with chromium sources will help understand the system- 
atic uncertainty and the detector cross section, reducing the uncertainty. It is also important 
to continue the gallium experiments to the statistics limit to establish consistency or incon- 
sistency with the 100 SNU benchmark. Theoretically, those models which predict 100 SNU 
can be compared with helioseismology data. In fact, some of the nonstandard models (the 
low T c model J37j, the large Sri model @, the mixing model |36[| , and the low Y model 



are in conflict with the sound speed profile inferred from helioseismology observations 
and therefore excluded. Further detailed testing of those nonstandard solar models with 
helioseismology data would be welcome. 

The flux constraints when the gallium and Kamiokande results are separately ignored are 
shown in Figs. ^| and |7|. The flux constraints from various combinations of the existing data 
are summarized in Table [TTT] and |TV| . We note that any combination of two experiments are 
consistent with the complete absence of the 7 Be and CNO fluxes, the 8 B flux of about 40% of 
the SSM, and the maximum pp flux, contradicting astrophysical solutions in general. That 
is, we have to ignore two of the three data to find a reasonable astrophysical explanation of 
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the solar neutrino problem. 



B. Future 

Since the current results are almost limited by systematic uncertainties, the present sta- 
tus described in the previous section is unlikely to change with the existing experiments 
unless there is a drastic revision in the data analyses. We expect, on the other hand, that 
our understanding of solar neutrinos will greatly improve once the results from the new gen- 



eration of high- statistic experiments are available. SNO [17| and Super-Kamiokande |18| will 
start in 1996, measuring the 8 B flux with high precision. The neutral current (NC) measure- 
ment in SNO and the measurements of the energy spectrum and time dependence in the two 
experiments will either confirm or rule out the neutrino oscillation hypothesis. BOREXINO 
f2(| will operate later in the decade and measure the 7 Be line spectrum separately. HELLAZ 
Plfl and HERON [22 can measure the pp and 7 Be neutrinos separately. 



Assuming that neutrino physics effects are absent, we should be able to calibrate solar 
models with precision measurements of the 7 Be and 8 B fluxes, independent of the existing 
experiments. The relevant flux parameter space with various standard and nonstandard 
solar models is displayed in Fig. ||(a). The determination of the 7 Be and 8 B fluxes at better 
than the 20% level should distinguish between competing solar models. It is also important 
to compare the future results to the present constraints for a consistency check among the 
data. 

For experiments sensitive to the pp flux, such as the gallium experiments, HELLAZ, 
and HERON, the relevant flux parameter space will be in Fig. |8|(b) and (c). The pp flux 
is, however, strongly constrained by the solar luminosity, and, to further distinguish the 
competing solar models, measurement uncertainties as small as a few % will be required. 

To study the sensitivity for determining the fluxes and distinguishing solar models, we 
have carried out a joint analysis assuming various possible outcomes from the new gen- 
eration experiments. The constraint on the fluxes is shown in Fig. ||(a) when SNO or 
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Super-Kamiokande data are assumed to be 0.50 ± 0.05 SSM. The allowed parameter space 
is for 90% C.L. The constraint on the 7 Be flux from BOREXINO Qdata (1.0 ±0.1 SSM) is 
displayed in Fig. |9](b). Displayed in Fig. [irj(a) are the constraints on both fluxes when the 
7 Be flux is measured in BOREXINO at the SSM value with an experimental uncertainty of 
10%; various values (0.3, 0.5, 0.7, and 1.0 SSM) for the 8 B flux measurement are assumed 
with a 10% experimental uncertainty. Fig. |To|(b) is the same except that various 7 Be values 
are assumed for a fixed central value of the 8 B flux (0.50 ± 0.05 SSM). 



The constraints are shown in Fig. |TT| when different experimental uncertainties are used 
for the 8 B flux measurements and for BOREXINO. With measurement uncertainties at 
the 10% level in SNO, Super-Kamiokande, and BOREXINO, the 7 Be and 8 B fluxes are 
determined accurately enough that the observations can distinguish between standard and 
nonstandard solar models and perhaps even constrain the SSM parameters. The constraints 
from the future data should be compared with the current constraint (Fig. 0) for a consis- 
tency check. 



III. FLUX CONSTRAINTS ASSUMING MSW 
A. Present 

Once the MSW effect is introduced in the analysis, the calibration of the neutrino fluxes 
becomes more complicated. One must constrain the initial fluxes and the MSW parame- 
ters simultaneously, while the neutrino spectrum can be distorted depending on the MSW 
parameters. We consider the simplest scenario, the two-flavor MSW effect. With the three 
existing data points and using Am 2 , sin 2 29, and 0(B) as completely free parameters, one 
obtains 

0(B)MB) SSM = 1-15 ± 0.53 (la), (3) 



The same analysis applies for the HELLAZ and HERON experiments. 
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while the other fluxes are fixed to the SSM values. Although the constraint is weak, it 
is consistent with the SSM predictions and already excludes (in the MSW context) some 
of the nonstandard models with a smaller 8 B flux. Since half of the SSM 8 B flux is seen 
in Kamiokande and since the MSW effect only reduces observed rates, the 8 B flux cannot 
be too small. Taking into account the Homestake and gallium data and also the neutral 
current contribution in Kamiokande, the 90% lower and upper limit is 0.47 and 2.07 of SSM, 



respectively. The constraint at 90% C.L. is displayed in Fig. [12[ The x distribution and 



the corresponding constraints on the MSW parameters are shown in Fig. £L3. 

If the 7 Be flux is introduced as an additional free parameter, no realistic constraint is 
obtained, even though the chlorine and gallium experiments have sensitivity to the flux. This 
is because the MSW survival probability for the 7 Be flux can be zero, allowing essentially 
any amplitude for the initial flux. In principle, the 7 Be flux has an upper limit due to the 
luminosity constraint, but the constraint is weak and irrelevant. We have repeated the fit 
by assuming smaller uncertainties for all experiments and by incorporating the Kamiokande 
spectral and day-night data |3]], but no constraint was obtained for the 7 Be flux. If MSW is 
operative, one needs a neutral current sensitivity for the flux (as in BOREXINO, HELLAZ, 
and HERON) to extract the 7 Be amplitude, which we will discuss in the next section. 

The core temperature, although model dependent, can be determined from the existing 
data using the power law for the 7 Be and 8 B fluxes. The power law obtained from a Monte 
Carlo estimation is 

0(Be) ~ T 8 and 0(B) ~ T™, (4) 

and the Tq dependence of the pp flux is obtained from the above relation and the luminosity 
constraint (Eq. [D), assuming the exponents of the pep and CNO fluxes as 2.8 and 22, respec- 
tively [[nj. The flux uncertainties from the nuclear reaction cross sections are included for 



S17 and S*34 as described in [I~J[]. The detector cross section uncertainties for chlorine and 
gallium are also included. As a result of a three parameter fit (two MSW parameters and 



Tc) [^g,|I4|]? we obtain 
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T c = 1.00 ± 0.03 (la) (5) 

in units of the SSM prediction (Tc = 1 = 15.57 x 10 6 K). The result is consistent with the 
SSM (Tc = 1 ± 0.006). That is, allowing the MSW effect, the present data determine Tc 
to within 3% and are consistent with the SSM predictions. We note that without the MSW 
effect no temperature could describe the data simultaneously [pREflBH] . The x 2 distribution 



and the constraints on the MSW parameters are shown in Fig. [14. 



B. Future 

In the next decade or so, the new generation of solar neutrino experiments will start and 
provide high-statistics data. Those experiments will measure the fluxes precisely and will 
allow a separation of the 8 B, 7 Be, and pp fluxes. Then, to determine the initial neutrino 
spectrum in the presence of nonstandard particle physics effects, what needs to be measured, 
and with what accuracy? 

We answer these questions quantitatively in the model independent framework, assuming 
two-flavor MSW oscillations, since it is the simplest solution of the solar neutrino problem 
and most successful in describing the existing data. We assume that the measurement 
of the charged to neutral current ratio in SNO will establish neutrino oscillations. We also 
assume that the measurement of the energy spectrum distortions and the day-night effect in 
SNO and Super-Kamiokande will distinguish the three separate MSW parameter branches 
from each other Jl4| and from vacuum oscillations |[39|| . The adiabatic and nonadiabatic 
regions will show 8 B spectrum depletion at higher and lower energies, respectively, which 
will be observable in SNO and Super-Kamiokande ||14|| . Most of the large-angle region shows 



the Earth effect, which will be measurable as day-night asymmetries or diurnal signal vari- 
ations in SNO, Super-Kamiokande, and BOREXINO OHO]. Since the spectrum and time- 



12 Similar analyses should be applicable to other particle physics effects if they do not involve too 
many new parameters. 
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variation information constrain the MSW parameters independent of the flux uncertainties, 
it would be best to incorporate those data directly in the analysis once the actual data are 
available. At present, however, we do not attempt to consider such constraints. We only 
consider the averaged SNO NC and BOREXINO rate (and the averaged Super-Kamiokande 
rate for some cases) as the minimal hypothetical data from the future experiments. Even 
so, one should be able to determine all of the parameters reasonably well. 

In the MSW calculations, we employ the electron density profile function and the neu- 
trino production profile functions of the Bahcall-Pinsonneault SSM. These functions are 
solar model dependent and should, in principle, be an additional source of uncertainties in 
constraining the fluxes and the MSW parameters. We have previously investigated those 
uncertainties by using three different SSMs and also by changing the peak location of the 
production profiles and the electron density scale height by 10% each [p4| . The effect on the 
obtained MSW parameters was negligible in the combined fit. 



In Figures |i5H19|, we consider the flux constraints for various possible outcomes of the 
SNO NC, BOREXINO, and Super-Kamiokande experiments that are consistent with the 
assumption that the MSW parameters are in the nonadiabatic (diagonal) branch. We in- 
clude the current results of the Homestake, (time-averaged) Kamiokande, and the combined 
gallium experiments, incorporating the detector cross section uncertainties in the radio- 
chemical experiments. However, omitting either the Homestake or gallium results does not 
change the result significantly, which will allow us to check consistency among data in the 
future. The constraints are obtained by fits to five free parameters [<f){pp), 0(Be), 0(B), 
Am 2 , and sin 2 29] imposing the luminosity constraint. The CNO and the minor fluxes are 
fixed to the SSM values. As shown later, using the CNO flux as an additional free parameter 
does not change our results significantly. 

Fig. |l5|(a) displays the constraints on the 7 Be and 8 B fluxes at 90% C.L. when the result 
of the SNO NC measurement is assumed to be the SSM value. The current data from Home- 
stake, Kamiokande, SAGE, and GALLEX are also included. The measurement uncertainties 
are taken as 10% of the signal. The SNO NC rate is unaffected by flavor oscillations, and 



yields a direct measurement of the 8 B flux. The 7 Be flux is not constrained even though 
the chlorine and gallium detectors have sensitivity, because the v e survival probability for 
this energy range can be zero, and therefore the initial 7 Be flux can take essentially any 
value. 

When the 7 Be measurement from BOREXINO Q is assumed to be 0.24 ± 0.024 of the 
SSM value, the allowed region is shown in Fig. p~J(b); the existing data are also included, 
but not the SNO result. Interestingly, both the 7 Be and 8 B fluxes are constrained in this 
case. The crucial factor is that, in electron scattering experiments, there is a contribution 
of the neutral current events from the or v T into which v e oscillates. For the 7 Be line 
spectrum at 0.862 MeV, the — e (or v T — e) cross section is 21% of the v e — e. Therefore 
there should be a signal of at least 21% of the initial flux even if the v e survival probability is 
zero. The BOREXINO result, combined with the existing data, gives a stringent constraint 
on both the original 7 Be flux and the MSW parameters, and those, in turn, constrain the 
8 B flux when combined with the Kamiokande result. 

The results when both the SNO NC and BOREXINO data are assumed for various 
different central values are shown in Fig. [TT5]. The MSW regions for the same or similar SNO 
NC and BOREXINO results are displayed in Fig. |Tj]. Considering that the constraints are 
independent of solar models, the allowed regions are determined surprisingly well. We also 
note that the information of the spectral distortion and of the day-night asymmetry, which 
is ignored here, will distinguish the adiabatic, nonadiabatic, and large-angle regions, and 



13 The charged current (CC) measurement in SNO, combined with the NC result, will determine 
the survival probability of the 8 B flux. Once neutrino oscillations are established, however, this 
information will not significantly improve the Kamiokande result included here. The effect of the 
CC measurement uncertainties are similar to the effect of the Super-Kamiokande measurement 
uncertainties discussed below. 

14 Our results apply for other 7 Be measurements with electron scattering, such as in HELLAZ 
and HERON. 
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therefore further constrain the MSW parameters. 

The effect on the flux constraints for various measurement uncertainties for SNO NC 



and BOREXINO are shown in Fig. |18[ Shown in Fig. |1!| is the constraint when the result 
from Super-Kamiokande is included for different measurement uncertainties. 

The effect of using the CNO flux as an additional free parameter is displayed in Fig. |2TJ; 
the result is essentially unchanged. 

The analysis has been repeated assuming measurements consistent with the MSW large- 
angle branch. The constraints from each SNO NC and BOREXINO measurement and the 



combined SNO NC and BOREXINO are displayed in Fig. ^T] and Fig. ^2|. The other results 
are essentially the same as for the nonadiabatic branch. 

From the analysis above, we conclude that, if the SNO NC and BOREXINO uncertain- 
ties are at the 10% level relative to signal, the 7 Be and 8 B fluxes should be constrained at 
the ± 20% and ± 15% level, respectively. This will clearly distinguish the standard and 
nonstandard solar models and perhaps even constrain the SSM parameters. The neutral 
current reaction for the 7 Be measurement ensures a non-zero signal (assuming flavor oscil- 
lations), which is especially important for obtaining stringent constraints on the neutrino 
fluxes and for distinguishing between competing solar models. 

Although model dependent, we have also carried out a simultaneous fit of Tc and Syj with 
hypothetical outcome from SNO NC and BOREXINO. The current data are also included. 
The constraint on T c and Sn is shown in Fig. [23]. T c and Syj will be simultaneously 
determined at the ±4% and ±20% level (90% C.L.) 

IV. CONCLUSION 

We have demonstrated that a model independent analysis using the four relevant fluxes 
(pp, 7 Be, 8 B, and CNO) as free parameters subject to the luminosity constraint is a feasible 
scheme for neutrino spectroscopy, and therefore for testing solar models. The analysis is 
viable with both standard and nonstandard neutrinos. 
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Assuming standard neutrinos, the existing experiments give a poor fit and essentially 
exclude any solar models. Even allowing this poor fit, there is no reasonable explanation 
for the following constraints from the data: the 7 Be and CNO fluxes are zero, and the 8 B 
flux is about 40% of the SSM prediction; the pp flux is the maximum value allowed by the 
luminosity constraint: 

0(pp)/0(pp)ssM = L089 " i- 095 ( 6 ) 
0(Be)/0(Be) SS M < 0.07 (7) 
0(B)/0(B) SSM = 0.41 ±0.04 (8) 
0(CNO)/0(CNO) SS m < 0.26, (9) 

where the uncertainties are at la; the CNO flux includes the 13 N and 15 neutrinos, which 
are varied with the same scale factor. When the constraints are expressed as absolute fluxes, 
one obtains 

(j){pp) = (6.53 - 6.57) x 10 10 cm^sec" 1 (10) 
0(Be) < 0.34 x 10 9 cm^sec" 1 (11) 
0(B) = (2.33 ± 0.23) x 10 6 cm^sec" 1 (12) 
0(N) < 1.28 x 10 s cm^sec" 1 (13) 
0(0) < 1.11 x 10 s cm _2 sec _1 . (14) 

This severe suppression of the 7 Be flux relative to the 8 B flux is inconsistent with any of 
the explicit nonstandard solar models. This problem is made even worse if is lower than 
the values usually assumed. Even discarding any one of the three data, the constraints are 
consistent with the above. 

When the two-flavor MSW effect is introduced in the analysis, the flux constraint from 
the current data is weak, but consistent with the SSM, sufficient to exclude the nonstandard 
models with too-small 8 B fluxes: 

0(B)/0(B) SS M = 1-15 ±0.53 (la) (15) 
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or 

0(B) = (6.54 ± 3.02) x 10 6 cm- 2 sec _1 (la) (16) 

No meaningful constraint is obtained if the other fluxes are introduced as free parameters. 

We have also considered the flux constraints in the presence of two-flavor MSW by assum- 
ing various outcomes from the next generation high-counting experiments. Of course, one 
can always consider more complicated particle physics effects, such as three-flavor oscilla- 
tions involving sterile neutrinos. Here, however, we consider the simplest scenario, expecting 
that, should two-flavor MSW be the case, it will be established as the most likely solution 
by the NC measurement in SNO, and by spectral distortions and day-night asymmetry mea- 
surements in SNO and Super-Kamiokande. Assuming hypothetical outcomes from the SNO 
NC and BOREXINO measurements with realistic uncertainties for this simplest scenario, 
we found that the 7 Be and 8 B fluxes will be determined at the ±20% and ±15% levels, mak- 
ing competing solar models distinguishable even if the MSW effect is operative. The MSW 
parameters will also be determined with sufficient accuracy independent of solar models. We 
emphasize that the neutral current sensitivity for 7 Be neutrinos in BOREXINO, HELLAZ, 
and HERON is essential for obtaining such constraints. We did not incorporate the infor- 
mation from the SNO CC rate, Super-Kamiokande rate, spectral distortions, or day-night 
asymmetry; those data should provide more stringent constraints on the MSW parameters 
as well as on the fluxes. The pp flux can be measured by the HELLAZ and HERON exper- 
iments, but a measurement uncertainty at the few % level is required to determine the flux 
more accurately than the luminosity constraint. 
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TABLES 



TABLE I. The standard solar model predictions of Bahcall and Pinsonneault (BP SSM) [|| and 
of Turck-Chieze and Lopes (TL SSM) along with the results of the solar neutrino experiments. 





BP SSM 


TL SSM 


Experiments 


Kamiokande a (10 6 cm 2 sec 1 ) 


5.69 ±0.82 


4.4 ±1.1 


2.89 ± 0.41 (0.51 ± 0.07 BP SSM) 


Homestake b (SNU) 


8 ± 1 


6.4 ± 1.4 


2.55 ± 0.25 (0.32 ± 0.03 BP SSM) 


SAGE c & GALLEX d (SNU) 


131.5 ±1 


122.5 ± 7 


77 ± 9 (0.59 ± 0.07 BP SSM) 



a The combined result of Kamiokande II and III (total of 1670 days) is 2.89 +0.22/-0.21 (stat) ± 
0.35 (sys) xlO 6 cm^sec" 1 §. 

b The result through June, 1992 (Run 18 - 124) is 2.55 ± 0.17 (stat) ± 0.18 (sys) SNU @. 

c The combined result of SAGE I and II (through January, 1993) is 74 +13/-12 (stat) +5/-7 (sys) 

SNU §. 

d The combined result of GALLEX I and II (30 runs, through October, 1993) is 79 ± 10 (stat) ± 
6 (sys) SNU §. 



TABLE II. To simplify the notation, we use the following neutrino fluxes as units. These 
reference fluxes correspond to the Bahcall-Pinsonneault standard solar model with the helium 
diffusion effect M. 



cm 2 sec 1 


<£(pp)ssm 


6.00 x 10 10 


4>(pep)ssM 


1.43 x 10 8 


4>(hep)ssM 


1.23 x 10 3 


0(Be)ssM 


4.89 x 10 9 


0(B)ssm 


5.69 x 10 6 


0(N)ssm 


4.92 x 10 8 


0(O)sSM 


4.26 x 10 8 


0(F)ssm 


5.39 x 10 6 
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TABLE III. The constraints on fluxes from various combinations of the current data with and 
without the MSW effect. The uncertainties are at la, and the fluxes are in units of the reference 
values defined in Table yj. The constraints are converted to absolute fluxes in Table IV. The upper 
limit on the pp flux (1.095) is due to the luminosity constraint. Without the MSW effect, we note 
that the constraints are consistent with each other even if any one of the three data is ignored, 
but are inconsistent with the SSM and nonstandard solar models, which generally suppress the 8 B 
flux more than the 7 Be flux. When the MSW effect is present, a reasonable constraint is obtained 
only for the 8 B flux. The obtained flux is consistent with the SSM prediction, albeit with a large 
uncertainty. 



PP 



7 Be 



<B 



CNO 



Constraints without the MSW effect 
Kam + CI + Ga 1.089 - 1.095 

Kam + CI 1.084 - 1.095 

Kam + Ga 1.085 - 1.095 

CI + Ga 1.082 - 1.095 



< 0.07 

< 0.13 

< 0.13 

< 0.16 



0.41 ±0.04 
0.42 ± 0.04 
0.50 ±0.07 
0.38 ± 0.05 



< 0.26 

< 0.38 

< 0.56 

< 0.72 



Constraints with the MSW effect 
Kam + CI + Ga < 1.095 



1.15 ±0.53 



TABLE IV. The same as Table III, but in units of absolute fluxes. The 13 N and 15 fluxes 



are varied with the same scale factor in the fits. 





pp a 


7 Be b 


8 B c 


13 N and 15 d 


Constraints without the MSW effect 










Kam + CI + Ga 


6.53 - 6.57 


< 0.34 


2.33 ± 0.23 


< 1.28, 


< 1.11 


Kam + CI 


6.50 - 6.57 


< 0.64 


2.39 ± 0.23 


< 1.87, 


< 1.62 


Kam + Ga 


6.51 - 6.57 


< 0.64 


2.85 ±0.40 


< 2.76, 


< 2.39 


CI + Ga 


6.49 - 6.57 


< 0.78 


2.16 ±0.28 


< 3.54, 


< 3.07 


Constraints with the MSW effect 










Kam + CI + Ga 


< 6.57 




6.54 ±3.02 







a In units of 10 10 cm~ 2 sec _1 . 
b In units of 10 9 cm _2 sec _1 . 
c In units of 10 6 cm _2 sec _1 . 
d In units of 10 8 cm _2 sec _1 . 
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FIGURES 

FIG. 1. The constraints on the 7 Be and 8 B fluxes when the Kamiokande, Homestake, and 
the combined gallium results are fit separately. For each point in this plane, the data are fit to 
the pp and CNO fluxes subject to the luminosity constraint. The fluxes allowed by the Homestake 
and gallium result are below the dotted and dashed line, respectively. The fluxes allowed by 
Kamiokande is between the dot-dashed lines. 

FIG. 2. The flux constraints obtained from the combined Kamiokande, Homestake, and gal- 
lium results. The constraints are shown for the (a) 7 Be- 8 B, (b) pp- 7 Be, and (c) pp- s B planes. 
The best fit parameters are 4>(pp) /4>(pp) ssm = 1-095, </>(Be)/(/>(Be)ssM = 0, 0(B)/(/>(B)ssm = 0.41, 

and 0(CNO)/0(CNO)ssm = (Table [nj), but this fit is P oor: Xmin/ 1 d - f - = 3 - 3 > which is ex- 
cluded at 93% C.L. Also displayed are the Bahcall-Pinsonneault SSM 90% region (BP-SSM) §, 



the Bahcall-Ulrich Monte Carlo SSMs @, the Turck-Chieze-Lopes (TL) SSM [|, and various 
nonstandard solar models (see the text). The observations are inconsistent with any of those stan- 
dard and nonstandard solar models. Smaller Sn values decrease only the 8 B flux [as indicated 
by the downward arrow in (a)], and aggravate the discrepancy between the combined data and 
nonstandard solar models. 

FIG. 3. The flux constraints from the Kamiokande and gallium data, but without the Home- 
stake result. The combined fit again indicates the larger suppression of the 7 Be flux relative to 
8 B, consistent with the constraint including the Homestake result (Fig. ||). The C.L. contours in 
(a) correspond to 0(CNO) = 0, while the nonstandard models within the 99% C.L. (the mixing 
models, the large model, and the small £34 model) predict non-zero CNO fluxes, aggravating 
the disagreement with the data. 

FIG. 4. The flux constraints from the Kamiokande and gallium results, but without the 
Homestake data, when the standard CNO flux is assumed. A non-zero CNO flux aggravates the 
disagreement between the data and solar model predictions. 

FIG. 5. The gallium experiment results, the SSM gallium rates, and the gallium rates of non- 
standard solar models which predict the 8 B flux consistent with or close to the 8 B flux observed 
in Kamiokande (see the text for details). The nonstandard solar models consistent with Kamio- 
kande predict the gallium rate Rq& ~ 100 SNU, contradicting the combined observed rate, 77 ± 9 
SNU. The MSW solution obtained from the combined Kamiokande and Homestake data predicts 
i? Ga < 100 SNU El, consistent with the data. 
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FIG. 6. The flux constraints from the Kamiokande and Homestake results, but without the 
gallium data. The constraints are consistent with those including the gallium data (Fig. (2j). 



FIG. 7. The flux constraints from the Homestake and gallium results, but without the Kam- 
iokande data. The constraints are consistent with those including the Kamiokande data (Fig. |2|). 

FIG. 8. Various standard and nonstandard solar models displayed in the (a) 7 Be- 8 B, (b) pp- 7 Be 
and (c) pp- s B flux parameter space. SNO and Super-Kamiokande will measure the 8 B flux. The 
SNO NC measurement will constrain the 8 B flux even if neutrino flavor oscillations are present. 
BOREXINO, HELLAZ, and HERON will measure the 7 Be flux. HELLAZ and HERON will also 
be capable of measuring the pp flux. The determinations of the initial 7 Be and 8 B fluxes at the 
< 20% level will make competing solar models distinguishable. For the pp flux, a determination at 
the few % level would be useful. 

FIG. 9. The flux constraints when the hypothetical results from (a) the 8 B flux measurement in 
SNO and Super-Kamiokande and (b) the 7 Be flux measurement in BOREXINO (and in HELLAZ 
and HERON) are considered. The standard neutrino properties are assumed. 

FIG. 10. The flux constraints for the combined SNO/Super-Kamiokande and BOREXINO 
results. The standard neutrino properties are assumed. The constraints are for (a) different 
SNO/Super-Kamiokande rates and (b) different BOREXINO rates. 

FIG. 11. The flux constraints for various measurement uncertainties in (a) 

SNO/Super-Kamiokande and (b) BOREXINO. The standard neutrino properties are assumed. 
With the measurement uncertainties at the 10% level, one can distinguish between standard and 
nonstandard solar models and perhaps even constrain the SSM parameters. 

FIG. 12. The flux constraint from the existing data when the MSW effect is assumed. The 
current data constrain 0(B) /0(B) SS m = 0.47 - 2.07 (90% C.L.) as shown in Fig. [13| (a). The solar 
models with too small 8 B fluxes are inconsistent with the existing data and the MSW hypothesis. 



The corresponding allowed MSW parameter space is displayed in Fig. 13 (b). No reasonable 



constraint is obtained when the fluxes other than 8 B are used as free parameters. 
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FIG. 13. The MSW-0(B) simultaneous fit to the existing data. This is a 3 parameter fit for 
8 data points, including 6 Kamiokande day-night data bins (5 d.f.) (a) The x 2 distribution as a 
function of 0(B). The current data constrain 0(B)/0(B) S sm = L15 ± °' 53 ( b ) The MSW 

allowed regions. The corresponding constraints on the 8 B flux are displayed in Fig. [12]. There is 
a third allowed region around sin 2 26 ~ 1 and Am 2 ~ 0.5 x 10~ 7 eV 2 , which is too small to be 
shown in the figure. Also shown is the region excluded by the Kamiokande day-night data (95% 
C.L., dotted line), which is independent of the 8 B flux uncertainty. For comparison, the allowed 
regions obtained assuming the Bahcall-Pinsonneault SSM and its uncertainties are also shown. 

FIG. 14. The MSW-Tc simultaneous fit to the existing data. This is a 3 parameter fit for 
8 data points, including 6 Kamiokande day-night data bins (5 d.f.) (a) The \ 2 distribution 
as a function of Tq- The current data constrain Tq = 1.00 ± 0.03, consistent with the SSM 
(Tq = 1± 0.006). (b) The MSW allowed regions. There is a third allowed region around sin 2 26 ~ 1 
and Am 2 ~ 0.7 x 10~ 7 eV 2 , which is too small to be shown in the figure. Also shown is the re- 
gion excluded by the Kamiokande day-night data (95% C.L., dotted line), which is independent of 
Tq- For comparison, the allowed regions obtained assuming the Bahcall-Pinsonneault SSM and its 
uncertainties are also shown. 

FIG. 15. The flux constraints for the MSW nonadiabatic region when the existing data plus 
possible results from (a) SNO and (b) BOREXINO are considered. The pp, Be, and 8 B fluxes are 
fit as free parameters subject to the luminosity constraint. These are 5 parameter fits (2 MSW 
parameters and 3 fluxes) to 4 data points (3 existing data plus 1 future data) with the luminosity 
constraint. 

FIG. 16. The flux constraints when the combined existing data plus possible results from both 
SNO and BOREXINO are considered. The projected experimental results are motivated by the 
MSW small-angle (nonadiabatic) solution. The constraints are for (a) different SNO NC rates and 
(b) different BOREXINO rates. These are 5 parameter fits (2 MSW parameters and 3 fluxes) to 
5 data points (3 existing data plus 2 future results) with the luminosity constraint. 

FIG. 17. The MSW allowed region when the pp, 7 Be, and 8 B fluxes are fit as free parameters 
with the luminosity constraint. The existing data plus the results from SNO NC and BOREXINO 
are used. We assume (a) different SNO NC rates with a fixed BOREXINO rate and (b) different 
BOREXINO rates with a fixed SNO NC rate. The constraints for the fluxes with similar assump- 



tions are shown in Fig. 16 and 22. Using the CNO flux as an additional free parameter does not 



change the allowed regions significantly. 
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FIG. 18. The flux constraints for the MSW nonadiabatic region for various measurement 
uncertainties in (a) SNO and (b) BOREXINO. The existing data are also included in the fits. 



FIG. 19. The flux constraint when hypothetical Super-Kamiokande results for various mea- 
surement uncertainties are included. The joint fit also includes the existing data and the hypo- 
thetical SNO NC and BOREXINO results. The MSW parameters are in the nonadiabatic region. 

FIG. 20. The flux constraints when the CNO flux is used as an additional free parameter. 
The joint fit includes the existing data and the hypothetical SNO NC and BOREXINO results. 
The MSW parameters are in the nonadiabatic region. 



FIG. 21. Same as Fig. 15, but the MSW parameters are in the large-angle region. 



FIG. 22. Same as Fig. 16, but the MSW parameters are in the large-angle region. 



FIG. 23. The constraints for Tc and Svj in the presence of MSW oscillations when the 
existing data plus both hypothetical SNO NC and BOREXINO results are considered. This is a 
4 parameter fit (2 MSW parameters plus Tc and Sn) to 5 data points (3 existing data plus 2 
hypothetical results) with the luminosity constraint. 
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FIG. 1. The constraints on the r Be and 8 B fluxes when the Kamiokande, Homestake, and the combined gallium results are fit separately. For each point in this plane, the 
data are fit to the pp and CNO fluxes subject to the luminosity constraint. The fluxes allowed by the Homestake and gallium result are below the dotted and dashed line, 
respectively. The fluxes allowed by Kamiokande is between the dot-dashed lines. 



33 



1.0 



CO 
CO 



PQ 



PQ 



0.5 



Monte Carlo SSMs 

• TL SSM 
OLow Z 

□ Low opacity 
■ WIMP 
A Large S n 
O Small S 34 
V Large S 33 

♦ Mixing models 

A Dar-Shaviv model 



V 




No oscillations 
Combined fit of 
current data 
90, 95, 99% C.L. 



BP SSM- 
♦ 



A 



Smaller S , 



T power law 



0.0 



0.0 



0.5 

4»( 7 Be) / 4»( ? Be) 



1.0 



1.0 



s 

CO 
CO 

PQ 



0.5 



PQ 



0.0 



Cjfc. ■ ^ BP SSM 


ii 
Monte Carlo SSMs 
• TL SSM 




OLowZ 


D * T c power law 

" + / 


□ Low opacity 
■ WIMP 
A Large S u 


A 


A Small S 34 
V Large S 33 


♦ ♦ 


♦ Mixing models 


A 


A Dar-Shaviv model 


v i 

V 

A 


.♦ 

\o 

A 

K ■ 


No oscillations 


\ \ 


Combined fit of 


\ K 


current data 




90, 95, 99% C.L. 



1.00 



SSM 



1.05 
4»(pp) / <|>(pp) 



1.10 



SSM 



(a) 7 Be- 8 B 



(b) pp- 7 Be 



34 



1.0 



CO 
CO 



PQ 



PQ 



0.5 



0.0 



"• * BPSSM 


Monte Carlo SSMs 
• TL SSM 

OLow Z 




□ Low opacity 
■ WIMP 
▲ Large S , , 

oil 


♦ 

□ 


O Small S 34 
V Large S 33 
♦ Mixing models 


♦ 

" A ▲ 


A Dar-Shaviv model 


T c power law — » . 

A V, 




▲ 

■ 


No oscillations 
\, * Combined fit 

90, 95, 99% C.L. 

▲ 



1.00 



1.05 
(|)(pp) / 4>(pp) SSM 



1.10 



(c) pp- 8 B 



FIG. 2. The flux constraints obtained from the combined Kamiokande, Homestake, and gallium results. The constraints are shown for the (a) 7 Be- 8 B, (b) pp- 7 Be, and (c) 
pp- s B planes. The best fit parameters are <f>(pp)/ '<f>(pp)ssM = 1.095, 0(Be)/0(Be) SS M = 0, <f>(B)/ <f>(B)ssM = 0.41, and 0(CNO)/0(CNO)ssm = (Table III), but this fit is poor: 
Xmin/1 d - f - = 3 - 3 > which is excluded at 93% C.L. Also displayed are the Bahcall-Pinsonneault SSM 90% region (BP-SSM) [8], the Bahcall-Ulrich Monte Carlo SSMs [23], the 
Turck-Chieze-Lopes (TL) SSM [9], and various nonstandard solar models (see the text). The observations are inconsistent with any of those standard and nonstandard solar 
models. Smaller Sir values decrease only the 8 B flux [as indicated by the downward arrow in (a)], and aggravate the discrepancy between the combined data and nonstandard 
solar models. 
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FIG. 3. The flux constraints from the Kamiokande and gallium data, but without the Homestake result. The combined fit again indicates the larger suppression of the 7 Be 
flux relative to 8 B, consistent with the constraint including the Homestake result (Fig. 2). The C.L. contours in (a) correspond to c/>(CNO) = 0, while the nonstandard models 
within the 99% C.L. (the mixing models, the large S33 model, and the small S34 model) predict non-zero CNO fluxes, aggravating the disagreement with the data. 
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FIG. 4. The flux constraints from the Kamiokande and gallium results, but without the Homestake data, when the standard CNO flux is assumed. A non-zero CNO flux 
aggravates the disagreement between the data and solar model predictions. 
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FIG. 5. The gallium experiment results, the SSM gallium rates, and the gallium rates of nonstandard solar models which predict the 8 B flux consistent with or close to the 
8 B flux observed in Kamiokande (see the text for details). The nonstandard solar models consistent with Kamiokande predict the gallium rate i?Ga ^ 100 SNU, contradicting 
the combined observed rate, 77 ± 9 SNU. The MSW solution obtained from the combined Kamiokande and Homestake data predicts i?Ga < 100 SNU [14], consistent with the 
data. 
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FIG. 6. The flux constraints from the Kamiokande and Homestake results, but without the gallium data. The constraints are consistent with those including the gallium 
data (Fig. 2). 
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FIG. 7. The flux constraints from the Homestake and gallium results, but without the Kamiokande data. The constraints are consistent with those including the Kamiokande 
data (Fig. 2). 
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FIG. 8. Various standard and nonstandard solar models displayed in the (a) r Be- 8 B, (b) pp- r Be and (c) pp- 8 B flux parameter space. SNO and Super-Kamiokande will 
measure the 8 B flux. The SNO NC measurement will constrain the 8 B flux even if neutrino flavor oscillations are present. BOREXINO, HELLAZ, and HERON will measure 
the 7 Be flux. HELLAZ and HERON will also be capable ol measuring the pp flux. The determinations ol the initial r Be and 8 B fluxes at the < 20% level will make competing 
solar models distinguishable. For the pp flux, a determination at the lew % level would be uselul. 
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FIG. 9. The flux constraints when the hypothetical results from (a) the 8 B flux measurement in SNO and Super-Kamiokande and (b) the r Be flux measurement in BOREXINO 
(and in HELLAZ and HERON) are considered. The standard neutrino properties are assumed. 
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FIG. 10. The flux constraints for the combined SNO/Super-Kamiokande and BOREXINO results. The standard neutrino properties are assumed. The constraints are 
(a) different SNO/Super-Kamiokande rates and (b) different BOREXINO rates. 
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FIG. 11. The flux constraints for various measurement uncertainties in (a) SNO/Super-Kamiokande and (b) BOREXINO. The standard neutrino properties are assumed. 
With the measurement uncertainties at the 10% level, one can distinguish between standard and nonstandard solar models and perhaps even constrain the SSM parameters. 
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FIG. 12. The flux constraint from the existing data when the MSW effect is assumed. The current data constrain 0(B)/c/>(B)ssm = 0.47 — 2.07 (90% C.L.) as shown in 
Fig. 13 (a). The solar models with too small 8 B fluxes are inconsistent with the existing data and the MSW hypothesis. The corresponding allowed MSW parameter space is 
displayed in Fig. 13 (b). No reasonable constraint is obtained when the fluxes other than 8 B are used as free parameters. 
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FIG. 13. The MSW-^(B) simultaneous fit to the existing data. This is a 3 parameter fit for 8 data points, including 6 Kamiokande day-night data bins (5 d.f.) (a) The \ 2 
distribution as a function of cf>(B). The current data constrain 0(B)/c/>(B)ssm = 1.15 ± 0.53 (lc). (b) The MSW allowed regions. The corresponding constraints on the 8 B flux 



are displayed in Fig. 12. There is a third allowed region around sin 29 ~ 1 and Am 



0.5 x 10 eV , which is too small to be shown in the figure. Also shown is the region 



excluded by the Kamiokande day-night data (95% C.L., dotted line), which is independent of the 8 B flux uncertainty. For comparison, the allowed regions obtained assuming 
the Bahcall-Pinsonneault SSM and its uncertainties are also shown. 
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(a) X 2 (T C ) (b) MSW regions 



FIG. 14. The MSW- Tc simultaneous fit to the existing data. This is a 3 parameter fit for 8 data points, including 6 Kamiokande day-night data bins (5 d.f.) (a) The \ 2 
distribution as a function of Tc- The current data constrain Tc = 1.00 ±0.03, consistent with the SSM (Tc = 1± 0.006). (b) The MSW allowed regions. There is a third allowed 
region around sin 2 28 ~ 1 and Am 2 ~ 0.7 x 10 -r eV 2 , which is too small to be shown in the figure. Also shown is the region excluded by the Kamiokande day-night data (95% 
C.L., dotted line), which is independent of Tc- For comparison, the allowed regions obtained assuming the Bahcall-Pinsonneault SSM and its uncertainties are also shown. 
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FIG. 15. The flux constraints for the MSW nonadiabatic region when the existing data plus possible results from (a) SNO and (b) BOREXINO are considered. The pp, 
and 8 B fluxes are fit as free parameters subject to the luminosity constraint. These are 5 parameter fits (2 MSW parameters and 3 fluxes) to 4 data points (3 existing data 
1 future data) with the luminosity constraint. 
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FIG. 16. The flux constraints when the combined existing data plus possible results from both SNO and BOREXINO are considered. The projected experimental results are 
motivated by the MSW small-angle (nonadiabatic) solution. The constraints are for (a) different SNO NC rates and (b) different BOREXINO rates. These are 5 parameter fits 
(2 MSW parameters and 3 fluxes) to 5 data points (3 existing data plus 2 future results) with the luminosity constraint. 
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FIG. 17. The MSW allowed region when the pp, r Be, and 8 B fluxes are fit as free parameters with the luminosity constraint. The existing data plus the results from SNO NC 
and BOREXINO are used. We assume (a) different SNO NC rates with a fixed BOREXINO rate and (b) different BOREXINO rates with a fixed SNO NC rate. The constraints 
for the fluxes with similar assumptions are shown in Fig. 16 and 22. Using the CNO flux as an additional free parameter does not change the allowed regions significantly. 
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FIG. 18. The flux constraints for the MSW nonadiabatic region for various measurement uncertainties in (a) SNO and (b) BOREXINO. The existing data are also included 
in the fits. 
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FIG. 19. The flux constraint when hypothetical Super-Kamiokande results for various measurement uncertainties are included. The joint fit also includes the existing data 
and the hypothetical SNO NC and BOREXINO results. The MSW parameters are in the nonadiabatic region. 



57 



2.0 



1.5 



MSW small-angle solution 
BOREXINO 0.24 + 0.024 SSM 
SNONC 1.0 + 0.1 SSM 4 

Standard CNO 

CNOfree 



T c power law 



S 

C/5 



PQ 



PQ 



1.0 




BP SSM 



0.5 



0.0 



0.0 



♦ 

♦ i 

ov / 

V ■ 

A 



▲ A 



90% C.L. 



0.5 1.0 
(K 7 Be) / 4»( ? Be) 



1.5 



2.0 



SSM 



FIG. 20. The flux constraints when the CNO flux is used as an additional free parameter. The joint fit includes the existing data and the hypothetical SNO NC and 
BOREXINO results. The MSW parameters are in the nonadiabatic region. 
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FIG. 21. Same as Fig. 15, but the MSW parameters are in the large-angle region. 
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FIG. 22. Same as Fig. 16, but the MSW parameters are in the large-angle region. 
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FIG. 23. The constraints for Tc and Sir in the presence of MSW oscillations when the existing data plus both hypothetical SNO NC and BOREXINO results are considered. 
This is a 4 parameter fit (2 MSW parameters plus Tc and Sir) to 5 data points (3 existing data plus 2 hypothetical results) with the luminosity constraint. 
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